arXiv:l509.00314v2 [quant-ph] 22 Jan 2016 


Quantum phase transitions detected by a local probe using time correlations and 

violations of Leggett-Garg inequalities 


F. J. Gomez-Ruiz^ J. J. Mendoza-Arenas^F. J. Rodriguez^, C. Tejedor^, and L. Quiroga^ 

^Departamento de Fisica, Universidad de los Andes, A.A. 4976, Bogota D. C., Colombia 
‘^Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom and 
^Departamento de Fisiea Teoriea de la Materia Condensada and Condensed Matter Physies Center (IFIMAC), 

Universidad Autonoma de Madrid, 28049, Spain 
(Dated: January 25, 2016) 

In the present paper we introduce a way of identifying quantum phase transitions of many- 
body systems by means of local time correlations and Leggett-Garg inequalities. This procedure 
allows to experimentally determine the quantum critical points not only of finite-order transitions 
but also those of infinite order, as the Kosterlitz-Thouless transition that is not always easy to 
detect with current methods. By means of simple analytical arguments for a general spin-1/2 
Hamiltonian, and matrix product simulations of one-dimensional XXZ and anisotropic XY models, 
we argue that hnite-order quantum phase transitions can be determined by singularities of the time 
correlations or their derivatives at criticality. The same features are exhibited by corresponding 
Leggett-Garg functions, which noticeably indicate violation of the Leggett-Garg inequalities for 
early times and all the Hamiltonian parameters considered. In addition, we hnd that the infinite- 
order transition of the XXZ model at the isotropic point can be revealed by the maximal violation 
of the Leggett-Garg inequalities. We thus show that quantum phase transitions can be identihed by 
purely local measurements, and that many-body systems constitute important candidates to observe 
experimentally the violation of Leggett-Garg inequalities. 


I. INTRODUCTION 


In recent years, quantum phase transitions (QPTs) of 
many-body systems have been the object of intense re- 
searctP^. This is the case not only due to the intrinsic 
interest that critical phenomena exhibit but also because 
the understanding and development of new states in con¬ 
densed matter or atomic systems may have prominent 
applications in areas such as high-temperature supercon- 
ductivit}!^ and quantum computatioiP. The seminal re¬ 
cent advances on quantum simulation schemed in sys¬ 
tems such as cold atoms in optical latticeJ^ and trapped 
ion^Ml constitute fundamental steps in this direction. 

Usually finite-order QPTs of a particular system are 
characterized by discontinuities of its ground state en¬ 
ergy or singularities of its derivatives with respect to the 
parameter that drives the transitions. Besides the deter¬ 
mination of order parameters, quantities such as gaps, 
spatial correlation functions, and structure factors are 
commonly used to determine the quantum critical points 
of several models. Remarkably, a few years ago it was 
realized that entanglement plays a fundamental role in 
critical phenomena and that different measures of entan¬ 
glement can be used to determine the location of several 
types of qptP^. Furthermore, the relation of Bell in¬ 
equalities and criticality has been recently explorecP^^^^. 

Since nonlocal measurements are not always accessible, 
in this paper we propose an alternative form to charac¬ 
terize QPTs by exploiting single-site protocols to obtain 
bulk properties of many-body systemJ^^Jt^. We argue 
that local time correlations can indicate the location of 
critical points for finite-order QPTs, in a similar way 
to measures of bipartite entanglement such as concur¬ 
rence and negativit^nH. This is exemplified by numeri¬ 


cal simulations, based on tensor network algorithms, of 
time correlations of one-dimensional (ID) spin-1/2 lat¬ 
tices described hy XXZ and anisotropic XY Hamiltoni¬ 
ans, which correspond to exhaustively-studied models of 
condensed matter physics. The first- and second-order 
transitions of these models are determined by nonana- 
lyticities of the time correlations and their first deriva¬ 
tive, respectively. We also relate QPTs to a different 
characterization of quantumness of a system, namely the 
violation of Leggett-Garg inequalities (LGI)P^^^, which 
indicates the absence of macroscopic realism and nonin- 
vasive measurability. We show that by maximizing the 
violation of these inequalities along all possible direc¬ 
tions, the infinite-order QPT of the XXZ model can be 
identified. Given that the models considered in our work 
describe several condensed-matter system^ and can be 
implemented in a variety of quantum simulator^, our 
analysis places them as interesting many-body scenar¬ 
ios for the experimental observation of the violation of 
Leggett-Garg inequalities. 


The paper is organized as follows: Section [TI| describes 
the ID spin models we focus on, whose QPTs are well 
known and thus allow us to check the adequacy of our 
proposal. Section [Hi] discusses how finite-order QPTs 
can be identified from local time correlations, provid¬ 
ing examples for the spin models previously described. 
Leggett-Garg inequalities and their role for determining 
QPTs are analyzed in Sec. El Finally, Sec. |V| contains 
the conclusions drawn from our work. 
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II. QUANTUM PHASE TRANSITIONS IN 
ONE-DIMENSIONAL SPIN-1/2 MODELS 

The main goal of our work corresponds to determin¬ 
ing whether local unequal-time correlations and Leggett- 
Garg inequalities can be used to localize and characterize 
QPTs. To provide an answer to this problem, we analyze 
the time correlations of systems described by spin-1/2 
time-independent Hamiltonians of the form 

^ = E E + E E (1) 

a. i,j a i 

Here erf denotes the Pauli operators at site i{a — z)^ 

is the coupling between spins at sites i and j along 
direction a, is the magnetic field at site i along direc¬ 
tion (a, and h = 1 . No restrictions on the dimensionality 
of the system or the range of the interactions are in prin¬ 
ciple required. 

While our analytical arguments are based on the 
Hamiltonian of Eq. 0 and thus are quite general (see 
Appendix , we restrict our numerical studies to two 
particular testbed Hamiltonians of condensed matter 
physics, namely the ID XXZ and anisotropic XY mod¬ 
els with nearest-neighbour interactions. These systems 
have been extensively studied in the literature, and their 
ground-state phase diagrams are very well knowiP^^^^. 
In this section we briefly describe the QPTs featured by 
these models. 


A. Spin-1 X X Z Model 

We first consider a ID system in which 1/2 spins are 
coupled through an anisotropic Heisenberg interaction. 
This case, known as the XXZ model, corresponds to 
4 ’^ = 4 ’^ = and = 0 . Thus it is 

described by the Hamiltonian 

H = JY^ [crfcrf_^i + (Tfcrf_^i + . (2) 

i 

Here the coupling J > 0 represents the exchange inter¬ 
action between nearest neighbors, and A is the dimen¬ 
sionless anisotropy along the 2; directioiP^. This model 
can be exactly solved by means of the Bethe ansat 
and possesses several symmetries. Namely, it features 
a continuous 1/(1) symmetry due to the conservation of 
the total magnetization in the 2; direction for any A and 
an additional SU{ 2 ) symmetry at A = ±1 due to the 
conservation of the total magnetization along the x and 
y directions. Furthermore, the Hamiltonian is invariant 
under transformations erf ^ —erf, thus having Z2 sym¬ 
metry. 

The model presents three different phases. First, for 
A < — 1 the ground state consists of a fully polarized 
configuration along the 2; direction, i.e. it corresponds 
to a ferromagnetic state. In the intermediate regime 
— 1 < A < 1 the system is in a gapless phase, which 


can be shown to correspond to a Luttinger liquid in 
the continuum limifP^. Finally, for A > 1, the ground 
state corresponds to an antiferromagnetic configuration. 
The ferromagnetic and gapless states are separated by a 
first-order QPT at A = —1, while the gapless and an¬ 
tiferromagnetic states are separated by a (infinite-order) 
Kosterlitz-Thouless QPT at A = 1. 


B. Spin-1 XY Model 


We now describe the anisotropic ID XY Hamiltonian 
for spins 1/2. It corresponds to = |J (1 -fy), = 

I J (1 — 7), Bl = Bz^ = B^ = By = 0 and is given by 


= jE 


l+7cc X , ~ ^ V V 




( 3 ) 

Here J > 0 represents the exchange interaction between 
nearest neighbors, 7 > 0 is the anisotropy parameter in 
the XY plane and > 0 is the magnetic field along 
the 2; direction. The limiting value 7 = 1 corresponds to 
the Ising model in a transverse magnetic field, which pos¬ 
sesses a Z2 symmetry, and the limit 7 = 0 is the isotropic 
XY model. In the thermodynamic limit A ^ 00, the 
anisotropic XY model can be exactly diagonalized by 
mea ns of Jordan-Wigner and Bogolyubov transforma- 
tion^SHMl. 


For the anisotropic case 0 < 7 < 1 the model belongs 
to the Ising universality class, and its phase diagram is 
determined by the ratio u = 2 BzlJ. When v > 1 the 
magnetic field dominates over the nearest-neighbor cou¬ 
pling, polarizing the spins along the 2; direction. This cor¬ 
responds to a paramagnetic state, with zero magnetiza¬ 
tion in the xy plane. On the other hand, when 0 < z/ < 1 
the ground state of the system corresponds to a ferromag¬ 
netic configuration with polarization along the xy plane. 
These phases are separated by a second-order QPT at 
the critical point u = 1 . Finally, for the isotropic case 
7 = 0 , a QPT is observed between gapless {u < 2 ) and 
ferromagnetic {u > 2 ) phases. We will only focus on 
the anisotropic model to illustrate the behavior of time 
correlations and LGI at criticality. 


III. SINGLE-SITE TWO-TIME CORRELATIONS 

Now we discuss how single-site two-time correlations 
(STG) can indicate different types of quantum phase 
transitions. We consider the symmetrized temporal cor¬ 
relation C{t) for a single-site operator A, given by 

C(t) = ^ {i’o\{A{t),A{ 0 )}\'ipo), ( 4 ) 

with 1 7/^0) fhe ground state of the time-independent 
Hamiltonian of interest iJ, the anticommutator be- 
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tween two operators, and A{t) the operator at time t, 

A{t) = ( 5 ) 

For simplicity, we consider that ^ 4 ( 0 ) corresponds to one 
of the Pauli operators of a particular site k [^(0) = , 

ji = z]. First note that the time correlations can be 

rewritten as 


Cit) = Re (^ol A( 0 )e-*^*^( 0 ) |^o)] , (6) 

with Eq the ground-state energy. Thus the correlation 
C{t) is a real quantity, which is fundamental for our sub¬ 
sequent analysis. As shown in Appendix for Hamilto¬ 
nian 0, the STC of Eq. 0 and their derivatives con¬ 
stitute appropriate quantities to determine the location 
of finite-order QPTs. In particular we obtain that the 
{p — l)th derivative of the STC is a function of Eq and 
its first p derivatives, so it can be written in the form 


dEo dPEo \ 
SAP-i V 9A ’■■■’ dXP A 


( 7 ) 


where A can be any Hamiltonian parameter, such as A for 
the XXZ model and u for the XY model. This means 
that, in general, a pth order QPT, which corresponds to 
a discontinuity or divergence of the pth derivative of the 
ground-state energy with respect to some Hamiltonian 
parameter, can be identified by the {p — l)th derivative 
of the STC with respect to the same parameter. Thus, 
a first-order QPT should in principle be identified by a 
discontinuity of the time correlations C(t), at any time 
Jt > 0, as a function of the parameter driving the transi¬ 
tion. Similarly, a second-order QPT should be recognized 
by a discontinuity or divergence of the first derivative of 
the time correlations with respect to the driving param¬ 
eter. Note that this result is similar to the observation 
of finite-order QPTs by measures of bipartite entangle- 
menfP^. Here, however, we are able to determine transi¬ 
tions by looking at a purely local (single-site) quantity. 

To provide stronger evidence that this is in fact the 
case, we calculate the time correlations for the XXZ and 
anisotropic XY models, and examine their behavior at 
the corresponding quantum critical points. Even though 
both models are exactly solvable, obtaining their physi¬ 
cal properties is a very challenging task. For example, for 
zero temperature exact time correlations are only known 
for A = 0 in the XXZ model (equivalent to the limit 
7 = 0 and = 0 in the XY model) and for A( 0 ) = a^SSl. 
Calculations based on a mean-field approach fail to repro¬ 
duce the time correlations correctly. Furthermore, exact 
diagonalization methods are restricted to small lattices. 
Thus to obtain quantitatively-correct results for much 
longer systems we perform numerical simulations based 
on tensor-network algorithms. Namely, we first obtain 
the ground state of both models for several parameters 
by means of the density matrix renormalization group 
algorithnPH, using a matrix product state descript iorP^. 
Subsequently we simulate the time evolution described 
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FIG. 1: STC along the z direction of the XXZ model, as a 
function Jt and anisotropy parameter A. The dashed green 
lines indicate the critical point A = — 1. 



A 


FIG. 2: STC along the z direction of the XXZ model, as a 
function of A, for Jt — 1 and Jt = 2. 


in Eq. ^ by means of the time evolving block decima- 
tiorP^. These methods allow us to carry out our sim¬ 
ulations efficiently, for lattices of several sites. In par¬ 
ticular, we consider systems of A = 100 spins (unless 
stated otherwise) with open boundary conditions, de¬ 
scribed by matrix product states with bond dimensions of 
up to X = 400 . Our implementation of the algorithms is 
based on the open-source Tensor Network Theory (TNT) 
librar}ffl. 


A. First-order QPT 

We start by observing the STC for the XXZ model, 
and focus on the transition between ferromagnetic and 
gapless states at A = — 1. All the results to be presented 
are in a time scale from 0 to 3 in units of 1/J. Since 
recent experiments on non-equilibrium spin models im¬ 
plemented in ultracold-atom quantum simulators have 
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been performed for similar {J/h = 27r x 8.6Hzf^ and 
even longer time scales {J/h = the effects we 

will present are in a time scale perfectly observable with 
current technology. In Fig. we show the correlations 
Cz{t) [i.e. for ^4(0) = cr|], evaluated at site k = 50, as a 
function of A and t. Additionally, in Fig. [^we plot the 
corresponding correlations at times Jt = 1 and Jt = 2 as 
a function of A. 

First, note that since the ground state of the system 
is ferromagnetic for A < —1, so IV^o) = IV^o) fhis 
regime, with the sign depending on the direction of polar¬ 
ization along the z axis. Furthermore, this state remains 
unchanged under magnetization-conserving time evolu¬ 
tion, such as that of the XXZ model. Thus the time 
correlations remain constant, with value Cz{t) = 1. For 
A > — 1 this is no longer the case. Since in this regime the 
states IV^o) are not fully polarized, they are strongly 
affected by time evolution. More importantly, when the 
system crosses the quantum critical point A = — 1 and 
enters the gapless state, the correlations exhibit a dis¬ 
continuous jump to values Cz{t) < 1 at any finite time 
Jt > 0, as depicted in Figs. and 

A similar result is obtained when calculating the STC 
Cx{t), i.e. with A(0) = cr^, which give identical results to 
the correlations along the y direction due to the symme¬ 
try of the Hamiltonian <§• These are shown in Fig. as a 
function of A and Jt, and in Fig. I^for two specific times, 
namely Jt = 1 and Jt = 2. In contrast to Cz{t), the cor¬ 
relations along the x direction do not remain constant 
in the ferromagnetic regime A < — 1, since flips the 
spin at site k and thus induces dynamics on the system. 
However, the correlations Cx{t) also show a discontinuity 
at A = — 1. Thus as expected from Eq. ([^, the different 
STC indicate the first-order QPT of the XXZ model by 
means of a discontinuity as a function of A at the quan¬ 
tum critical point. 

Note that neither Cz{t) nor Cx{t), or any of their 
derivatives, indicate the existence of the Kosterlitz- 
Thouless QPT at A = 1, given that it is of infinite order. 



FIG. 4: STC along the x direction of the XXZ model, as a 
function of A, for Jt — 1 and Jt — 2. 


However, we will observe in Section IV B that it is possi¬ 
ble to identify this transition by the maximization of the 
violation of LGI. 


B. Second-order QPT 

Now we consider the transition between ferromagnetic 
and paramagnetic phases of the anisotropic XY model. 
In particular, we illustrate the transition for two cases, 
namely the limit 7 = 1, which corresponds to the Ising 
model with a transverse magnetic field, and the inter¬ 
mediate case 7 = 0.5. We verified that the STC along 
any direction a = x^y^z give the same qualitative infor¬ 
mation regarding the QPT, which is in agreement with 
Eq. 0 - so we focus on Cz (t) and do not show the results 
of the other directions. 

In Eig. [^we show the ^ time correlations as a function 
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FIG. 3: STC along the x direction of the XXZ model, as a FIG. 5: STC along the z direction of the anisotropic XY 
function of Jt and anisotropy parameter A. model for 7 = 1 (Ising model), as a function of Jt and v. 
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FIG. 6 : STC (upper panel) and first derivative (lower panel) 
along the z direction of the anisotropic XY model (7 = 0.5,1), 
as a function of ly, for Jt — 1 and Jt — 2. 


of Jt and z/, for 7 = 1; the results for 7 = 0.5 are qualita¬ 
tively similar. In addition, we depict in the upper panel 
of Fig. 1^ the correlations for times Jt = 1 and Jt = 2 
as a function of the magnetic field, for both 7 = 1 and 
7 = 0.5. In contrast to the XXZ case, here the corre¬ 
lations are continuous for the whole range of values of 
u considered. However, the first derivative with respect 
to u is not a well-behaved function. As exemplified in 
the lower panel of Fig. [^for two particular times, 
shows a sharp maximum at the quantum critical point 

= 1. Thus, in accordance with Eq. Q, the second- 
order QPT of the model can be identified by means of a 
singularity in the first derivative of the local time corre¬ 
lations with respect to the Hamiltonian parameter which 
drives the transition, which in this case is u. 

IV. LEGGETT-GARG INEQUALITIES 

Since the birth of quantum mechanics, its non- 
deterministic nature and nonlocal structure have moti¬ 
vated many theoretical debates that have recently moved 
to the experimental field. In particular, Bell inequalities 
establish a natural border to the spatial quantum corre¬ 
lations in separate systems. Leggett and Gar^^ in 1985 



Jt 


FIG. 7: Leggett-Garg functions for the XXZ model. Upper 
panel: K^{t). The regions with diagonal lines corresponds to 
the regime of anisotropies A and Jt in which the x Leggett- 
Garg inequalities are violated. The white region below A = 
— 1 indicates that there the time correlations remain constant. 
Lower panel: K^{t). The regions with diagonal lines region 
corresponds to the regime of anisotropies A and Jt in which 
the X Leggett-Garg inequalities are violated. 


showed that the temporal correlations obey similar in¬ 
equalities. 

In our intuitive view of the world, probabilities are due 
to our uncertainty about the state of a system, but they 
are not a fundamental description of it. For example, 
when we toss a coin to the air, it has probability one 
half of landing tails or heads. We also assume that if we 
had the precise knowledge of its position and momentum, 
and enough computational power, we would be able to 
determine on which side the coin will land. We do not 
think that the coin is in a superposition of states, such a 
Schrodinger’s cat. This is known as macroscopic realism. 
In addition, we assume that making measurements on a 
system does not modify its present state, in the way pro¬ 
jective quantum measurements do. This is referred to as 
noninvasive measurability. Based on these two principles 
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FIG. 8: Maximum violation of Leggett-Garg inequalities for 
the XXZ model, along both z and x directions. 


Leggett and Garg obtained a set of inequalities, which 
is consistent with the macroscopic intuition. One form 
these LGI can take is 


C (ti, ts) - C (ti, t2) - C (t2, ts) > -1, (8) 

where C is the two-time cor¬ 

relation of a dichotomic observable Q (with eigenvalues 
g = ±1) between times ti and tj, and ti < t 2 < On 
the other hand, if the correlation functions C are 

stationary, i.e., they only depend on the time difference 
T = ti — tj, then the Leggett-Garg inequality (i can be 
written as^ 



II2.0 

--1.0 

--0.0 

—-1.0 



K_ (r) = C (2r) - 2C (r) > -1, (9) 

which defines the Leggett-Garg functions K_ (r) for time 
r. Just as with Bell inequalities, any system that vio¬ 
lates inequality shows some behavior that is essen¬ 
tially nonclassical. This is why violations of LGI are 
used as a measure of quantumnesJ^. In the following 
we discuss different Leggett-Garg functions cor¬ 

responding to measurements of spin components along 
the a direction, and see whether they can give informa¬ 
tion about the QPTs previously discussed. 


A. Finite-order QPTs 


We start by showing how the Leggett-Garg funct ions 
K^{t) signal the finite-order QPTs discussed in Sec.llII 


In Fig.we depict both K^{t) (upper panel) and K'_[t) 
(lower panel) for the XXZ model as a function of A and 
time. Regarding the results along the a = z direction, 
we first note that for A < —1 the value of the Leggett- 
Garg function remains equal to K^{t) = —1 for any time. 
Thus in the ferromagnetic phase the corresponding LGIs 
are never violated. This is clearly a direct consequence 


FIG. 9: Leggett-Garg functions K^{t) for the anisotropic 
XY model. Upper panel: 7 = 1 . Lower panel: 7 = 0.5. 
The regions with diagonal lines corresponds to the regime of 
parameter u and Jt in which the z Leggett-Garg inequalities 
are violated. 


of the constant value of the time correlations previously 
discussed (see Figs. and |^, and manifests the classi¬ 
cal nature of the ferromagnetic state when undisturbed. 
The situation is entirely different for A > —1. Not only 
K^(t) does vary on time, but it indicates a violation of 
the Leggett-Garg inequalities for early times. As to the 
results along the a = x direction, all the values of A 
considered show a violation of the inequalities for early 
times. In addition, the violation lasts longer as |A| de¬ 
creases. In Fig. we plot the maximum value we 

obtain for the violation of the Leggett-Garg inequalities 
as defined by 


= min KZ{t), (10) 

for both directions a = as a function of A. This 
clearly shows that similarly to time correlations, the first- 
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order QPT of the XXZ model can be identified by a 
discontinuity of the I/^ax function at the critical point. 
Also, the maximal violation occurs along the 2 : direction, 
close to the noninteracting limit A = 0. In contrast, for 
the magnetically-ordered phases, the maximal violation 
occurs along the x direction. 

The Leggett-Garg functions also help determine the 
second-order QPT of the anisotropic XY model. In 
Fig. [9 we show K^{t) for several values of as a func¬ 
tion^ time, for 7 = 1 (upper panel) and 7 = 0.5 (lower 
panel). Notably, for all the values of u considered, the 
system features the violation of the inequalities. Initially, 
for Jt < 0.5, the violation of the inequalities lasts longer 
as decreases. Interestingly, for longer times, revivals of 
the violations are seen for low values of u. Thus weak 
magnetic fields favor the observation of the violation of 
the Leggett-Garg inequalities along 2 ; direction. 

Just as the time correlations, the Leggett-Garg func¬ 
tions K^(t) and the maximal violation functions 
(see upper panel of Fig. 10) are continuous in the whole 
parameter regime. However, their first derivative tends 
to diverge at the quantum critical point z^ = 1 as the size 
of the system increases (see inset lower panel of Fig. 10). 
This is shown in the lower panel of Fig. 10 for and 
both 7 = 1 and 7 = 0.5. As expected, the behavior of the 
time correlations is translated to the Leggett-Garg func¬ 
tions, and they are able to signal the second-order QPT 
of the anisotropic XY model by means of a singularity 
in their first derivative. 



FIG. 10: Upper panel: Maximum violation of Leggett-Garg 
inequalities for the anisotropic XY model (7 = 1,0.5) along 
the 2 : direction, as a function of v. Lower panel: First deriva¬ 
tive with respect to h'. Inset: scaling of maximum of the 
derivative with respect to v 


B. Infinite-order QPT of the XXZ model 

We have observed that finite-order QPTs can in princi¬ 
ple be determined by means of a singular behavior of local 
unequal-time correlations and Leggett-Garg functions, or 
of their derivatives. However, this form is not suitable 
to identify infinite-order transitions. In fact, the results 
shown so far do not feature any singular property at the 
quantum critical point A = 1 of the Kosterlitz-Thouless 
transition of the one-dimensional XXZ model. However, 
it is possible to locate this transition from Leggett-Garg 
functions, as we discuss in the present Section. This is 
very similar to the observation of the transition from Bell 
inequalities^, with the notable difference that here we 
actually have violation of the respective inequalities, and 
thus we can perceive the quantumness of the system. 

The first point to note is that to actually establish 
that a violation of the inequalities exists, and also when 
the maximal violation occurs, we must consider all the 
possible directions a of evaluation of time correlations. 
For the XXZ model, this corresponds to a = z and 
a = X, the latter giving the same results as for a = y due 
to the symmetry of the Hamiltonian. For this we define 
the function 

-^max max (ff) 

a=z,x 


which maximizes over all times and directions the vio¬ 
lation of the Leggett-Garg inequalities. We show 
as a function of A in Fig. [m note that it indicates the 
first-order QPT at A = —1 by means of a discontinuity. 

The second point to note, responsible for the observa¬ 
tion of the Kosterlitz-Thouless transition by means of 
Bell inequalities, is that at the isotropic point of the 
Hamiltonian there is a change in the largest type of spa¬ 
tial correlatioiP^. Namely, for |A| > 1 and spins sepa¬ 
rated by r lattice sites, \ < |(crfcrf^^)|, while for 

-1 < A < 1 we have that \ > |(crfcrf^^)|. As ex¬ 

pected from the discussion of Appendix this behavior 
is translated to the local time correlations and related 
functions. In fact, as seen in Fig. ^max < imax 
the ordered phases, while for the gapless 

regime. As shown in Fig. mi a sharp local maximum ap¬ 
pears at A = 1 in the function, and a singularity of 
its first derivative results. Thus, by means of a function 
characterizing the total maximal violation of Leggett- 
Garg inequalities, we are able to locate the infinite- 
order Kosterlitz-Thouless transition of the XXZ model. 
It would be interesting to observe whether Kosterlitz- 
Thouless transitions for other quantum systems can be 
identified in this form. 













FIG. 11: Total maximum violation of Leggett-Garg inequal¬ 
ities for the XXZ model as a function of A. The light-red 
zones indicate the regimes in which the maximal violations 
comes from inequalities along the x direction, while the light- 
blue zone in between shows the regime in which the maximal 
violations occur along the z direction; see Fig. 


importantly, we found that at least for one direction, the 
Leggett-Garg inequalities are violated for early times and 
the whole regime of parameters considered, in contrast 
to Bell inequalities^. Furthermore, the maximization 
of this violation allowed us to identify the infinite-order 
Kosterlitz-Thouless QPT of the XXZ model at A = 1, 
which was not possible from the separate observation of 
time correlations along each direction. Given the large 
amount of materials desc ribed by the test-bed models 
discussed in our worlP^, and the seminal advances on 
their implementation in quantum simulator^, we expect 
that our results extend the range of systems in which the 
violation of Leggett-Garg inequalities can be observed 
experimentally^. 

For future research, it would be interesting to ana¬ 
lyze whether local time correlations and Leggett-Garg 
inequalities can identify the existence of different-order 
nonequilibrium quantum phase transitionJ^^^^. Similar 
analysis on Bell-Leggett-Garg inequalities^ could lead to 
important insights into the relation between measures of 
quantumness and quantum criticality. 


V. CONCLUSIONS 
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In the present paper we have discussed whether single¬ 
site time correlations and Leggett-Garg inequalities al¬ 
low the identification of QPTs in many-body quantum 
systems. By means of efficient matrix product simula¬ 
tions and analytical arguments, we have answered this 
question in the affirmative for different spin-1/2 models, 
for both finite- and infinite-order QPTs. Thus we have 
shown that QPTs can be detected by purely-local mea¬ 
surements. 

Initially, by means of a first-order approximation for a 
general spin-1/2 Hamiltonian, we argued that apth order 
QPT can be located by a singular behavior of the (p—l)th 
derivative of the local time correlations at the quantum 
critical point. Thus, these correlations indicate quan¬ 
tum criticality in a form similar to different measures of 
bipartite entanglemenff^. Furthermore, this behavior is 
directly transferred to the corresponding Leggett-Garg 
functions. 

To support this general result, we calculated several 
time correlations for large one-dimensional XXZ and 
anisotropic XY spin systems, using the density matrix 
renormalization group and time evolving block decima¬ 
tion methods. In particular, we showed that the first- 
order ferromagnetic-gapless QPT of the XXZ model 
is manifested as a discontinuity of the correlations at 
A = —1, along any possible direction and for any fi¬ 
nite time. Subsequently we showed that the second-order 
paramagnetic-ferromagnetic QPT of the anisotropic XY 
model is observed by means of a divergence of the first 
derivative of the correlations with respect to the magnetic 
field at = 1. 

We also showed that the Leggett-Garg functions can 
help identify finite-order QPTs in a similar fashion. More 
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Appendix A: Proof of the relations between 
two-time correlations and ground-state energy 


In this section we show how STC can indicate finite- 
order QPTs of spin-1/2 systems of any dimensionality or 
coupling range. Thus we consider a Hamiltonian of the 
form 




a i,j 


ri i T~>i 

=EE'^«"^+EE^o 


(Al) 


a 


‘dBi 


As shown in Eq. the time correlations of a single-site 
operator A are given by 

C{t) = Re (^ol A(0)e-'"M(0) |^o)] (A2) 

Now we expand the time evolution operator as 

oo / . xZ 

e-iHt ^ (A3) 

1=0 
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So the product ^4(0)6 *^^^4(0) can be written as 

OO /_ .,\Z 

^(0)e-*"M(0) = ^^A{0)H^A{0) 

1=0 

OO /_ -,\l 

= Y, ^-Jp- (^(O)ii'^(O))' = (A4) 

1=0 

where we have used 74(0)A(0) = T repeatedly in 

the first line to obtain the secon d equality (/i = x, z). 

Using the explicit form ( |A1| ) of the Hamiltonian, we 
obtain after straightforward algebra that 

H(0)H^H(0) = H- f^. (A5) 

where the second term explicitly depends on site k where 
the correlations are calculated, namely 


/fe = 2 E E (A6) 

\ j J 

Therefore, the STC is given by 

C{t) = Re (V’ol |V’o)] • (A7) 

To easily observe the relation between QPTs and time 
correlations, we restrict to first order in t he exponential 
within the expectation value of Eq. (A7). In this case, 
we have 


C{t) cos {Eot) + sin {Eot) [Eot - i (■*/'o| fk I'^/’o)] (AS) 
Now we rewrite in the form 




dJi 


,k 




\ 3 

Finally, using the Hellmann-Feynman relations 


V’o 


dH 


dX 


dE: 


(A9) 


^0 ) = ^, (AlO) 


with A any parameter of the HamiltoniarP^, we obtain 
that up to first order in the expansion of the time evolu¬ 
tion operator, the time correlations are given by 


identified by discontinuities of the STC as a function of 
Hamiltonian parameters. 


Now we consider the first derivative of Eq. (All) 
with respect to some Hamiltonian parameter, e.g. . 
We obtain 

dC{t) 




t cos (Eot) 


5Fo 


Eot-2tY 

\ 3 

- 2t sin {Eot) 

L j 

d^Eo ^ 


+ 4’" 







(A12) 


This means that the first derivative of the STC with re¬ 
spect to a Hamiltonian parameter is proportional to the 
second derivative of the ground-state energy with respect 
to the same parameter. 


dC{t) 


(X t 


d^Eo 




and 


dc{t) 


(X Eot. (A13) 


As a res ult of the first proportionality relation of 
Eq. (A13), the derivatives of unequal-time correlations 


indicate second-order QPTs by means of a discontinu¬ 
ity or divergence at the corresponding quantum critical 
points. 

The previous results indicate that, in general, a finite- 
order QPT can be identified by the properties of the STC 
of the system. Namely, given that 


and 


dp-^c{t) 

dP-^C{t) 


(X t 


(X t 


dPEo 

dP-'^Eo 


(A14) 


, dP-^C{t) ^ ^ 
d{Jp )P 1 


C{t) « cos {EqI) + sin {EqI) 



- 2 fE 



dEo 


+ B’l 


dEo\ 

dB^a ). 


(All) 


Thus we have obtained that the STC are proportional 
(apart from a structureless term Eot) to the first deriva¬ 
tives of the ground-state energy of the system, which 
show a discontinuity at the critical point of a first-order 
QPlEl. This means that first-order QPTs are directly 


the {p — l)th derivative of the STC with respect to 
some Hamiltonian parameter is a function of all gth 
derivatives of the ground state energy with respect to 
the same parameter, for 0 < g' < p. Thus a pth order 
QPT, which corresponds to a discontinuity or divergence 
of the pth derivative of the ground-state energy, can be 
identified by the (p — l)th derivative of the STC. 


The validity of this result is not affe cted by taking the 
expansion of the exponential of Eq. (A7) to higher or¬ 
ders. Consider, for instance, the second-order correction 
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which adds the terms 
C^^\t)=-^^cosiEot)x 


{E^, - 2Eo (tAol /r l^o) + (V^oK/D" l^o)) 

(A15) 


to the time correlations in Eq. (All). The term 


(V^ol fkli^o) has the form already displayed in Eq. (All). 


The third component (V^oK/^)^ IV^o) results in more com¬ 
plicated (up to three-site) expectation values in addition 
to more terms dEo/dJ^^ and dEoldBj^. These elements 
will continue appearing in higher-order expansions, either 
separately or in expectation values {^o\ fj^ \ipo)- So these 
expansions would lead to the observation of finite-order 
QPTs as previously discussed for the first-order case. 
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